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The two-fold singularity has played a significant role in our understanding of uniqueness 
and stability in piecewise smooth dynamical systems. When a vector field is discontinuous 
at some hypersurface, it can become tangent to that surface from one side or the other, and 
tangency from both sides creates a two-fold singularity. The local flow bears a superficial 
resemblance to so-called folded singularities in (smooth) slow-fast systems, which arise at 
the intersection of attractive and repelling branches of slow invariant manifolds, important 
in the local study of canards and mixed mode oscillations. Here we show that these two 
singularities are intimately related. When the discontinuity in a piecewise smooth system 
is regularized (smoothed out) at a two-fold singularity, the resulting system can be mapped 
onto a folded singularity. The result is not obvious since it requires the presence of nonlinear 
or ‘hidden’ terms at the discontinuity, which turn out to be necessary for structural stability 
of the regularization (or smoothing) of the discontinuity, and necessary for mapping to the 
folded singularity. 


I. INTRODUCTION 


If a flow is piecewise smooth, having a vector field that is discontinuous on some hypersurface 


S, then under generic conditions there can exist isolated singularities where the flow curves (or 
‘folds’ parabolically) towards or away from S on both sides of the surface. The result is a two-fold 
singularity, as depicted in figure [H generic in systems of three or more dimensions. 


Consider a piecewise smooth system of the form 



f(x; -1-1) if xi > 0 , 
f(x; —1) if xi < 0 , 


( 1 ) 


where f is differentiable with respect to the variable x = {xi,X 2 ,xs) A two-fold singularity is a 
point Xp = {xip,X 2 p,xsp) at which 


0 = xi = lim xx\^ / lim 


( 2 ) 


2 



FIG. 1: Three kinds of two-fold. The main figure shows the phase portrait: in the unshaded regions the flow 
crosses (cr.) through a discontinuity at xi = 0, in the shaded regions the flow can only slide along the discontinuity 
on Xi = 0, the region being attracting {a.si.) or repelling {r.sL). In the examples shown, determinacy-breaking 
occurs at the singularity, meaning that the flow there becomes set-valued, the set has 2 dimensions in (i) and 3 
dimensions in (ii-iii). 


and for which the two curves defined W lim^^o ®ilxi=±<5 = 0 on xi = 0 are transversal. A more 
complete description can be found in 

For a local singularity that is so easy to define, the ‘two-fold’ singularity has proven surprisingly 


difficult to characterize, from its first description in [a, [IT]] in three dimensions, to its study in higher 
dimensions in Q]. These exclusively consider the class of Filippov system obtained by expressing 
m as the convex combination 


X = f(x; A) = -hi) -b ^-^f(x; -1) , 


( 3 ) 


A G 


( 4 ) 


2 ' > ' ^ ' 2 

where 

sign (xi) if xi 7^ 0 , 

[—1, -1-1] if xi = 0 . 

The dynamics on the discontinuity surfaces in figure [T] arise, for example, from this expression. 
Early questions about the structural stability of two-fold singularities in such systems have been 
resolved by uncovering their intricate phase portraits, revealing various topologically stable phase 
portraits separated by bifurcations j^, 7, , which include the birth of limit cycles, bifurcation of an 

invariant nonsmooth diabolo, and passage of equilibria through the singularity. In its structurally 
stable forms (i.e. away from any bifurcations), the two-fold is neither an attractor nor repellor, so 
the flow either misses the singularity, or traverses it in finite time. 
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It is in the latter that important questions remain concerning the two-fold singularity, partic¬ 
ularly in cases where it forms a bridge from an attracting region on the discontinuity surface into 
a repelling region, mimicking canard behaviour of smooth two-timescale systems l|. In a loose 
definition suitable to both smooth and nonsmooth dynamics, canards are trajectories that persist 
from an attracting invariant manifold to a repelling invariant manifold. Numerous canards can be 
seen in figure [H Such behaviour takes a more extreme form in nonsmooth systems, because the 
two-fold breaks determinacy in both forward and backward time through the singularity. This is 
illustrated in figure [U where a typical single trajectory is shown entering the singularity, being 
deterministic until it does so, and afterwards exploding into a set-valued flow of infinite onward 
trajectories. This shape of this outset of the flow is determined by the local vector fields. 

Particularly because of some similarity to canard dynamics, attention has turned to how the 
two-fold can be understood as a limit or approximation of a smooth flow. An equivalence between 
“sliding” motion along a discontinuity surface, and “slow” motion on invariant manifolds of a 
smooth two-timescale system, has been shown |l4l.ll6l|. Concerning two-folds, the relation between 


the sliding phase portraits at two-: 


have received growing attention 


bid singiTarities, and canard dynamics of two-timescale systems. 


[3,4, 12 , 3 , 20 ]. In 0,4 


20 ( 1 . In , a qualitative association was made with 


the so-called ‘folded’ singularities of two timescale systems. 

A folded singularity can be defined in a system 

(£Xi,X 2 ,X 3 ) = (gi(xi, X 2 , X 3 ; e), g 2 (xi, X 2 , X 3 ; e), gsixi, X 2 , X 3 ; e)) , 

where the gi are differentiable, with 0 < e <C 1, as a point satisfying 

dgi 


(5) 


0 = gi = gi = 


dgi , , d^gi 


( 6 ) 


dxi ’ <9x2,3 ' " ' dx\ 

We shall prove here a more direct connection between the two singularities, by showing equivalence 
between the singularities defined by ([2]) and ([6]) under explicit coordinate transformations. 

This approach offers a different viewpoint on desingularizing the two-fold singularity to that 
taken in 112. 


1, 13, 3|- There, blow up methods are used to show that a regularization (a smoothing) 
of the two-fold contains canards. The results there apply specifically to the Sotomayor-Teixeira 
regularization [ 3 |- essentially replacing the sign function A in the convex combination ([3]), with a 
smooth sigmoid function of xi. While not currently in common use, it is easy to show that more 
general forms of system are possible, replacing ([3]) with the piecewise smooth systems 


X = f(x; A) = + 1 ) + ^^f(x; -1) -k (1 - A) , 


(7) 
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using Q, which like (l3|) coincides with ([T]) for xi 7 ^ 0 . The function g is an arbitrary finite vector 


3 , 0,121 


field, and the possibility of nonlinear A dependence can be found discussed already in 
(including an experimental model in {l3|). Sometimes called ‘hidden’ terms because they vanish 
everywhere except at the discontinuity, a general approach for handling nonlinear dependence on A 
was introduced conceptually in a. The S„t„.e.„,-Teixehe .heo.. of .....afion was e.fea.ed 
to m in {l3 |. This provides regularized systems that are not included in the Sotomayor-Teixeira- 
Filippov approach through ([3|), an issue we will illustrate with a simple example in section rVIIl 
Such systems turn out to be essential for the equivalence we seek to prove here. 

In section m we introduce the normal form of the two-fold singularity, and outline the basic 
steps for its study by regularizing the discontinuity in section IIIII In section IIVI we regularize 
the normal piecewise smooth system, assuming only linear dependence on A as in the standard 
literature, but show that this results in a degenerate system. In section |V] we perturb this using 
nonlinear dependence on A, finding that it breaks the degeneracyy, and can be mapped onto the 
folded singularity of a smooth two timescale system. Remarks showing that these results follow 
also if we blow up, rather than regularize, the discontinuity, are given in section IVlIl 


II. THE TWO-FOLD SINGULARITY 


The normal form of the two-fold singularity is 

j (-X2,ai,6i) if xi > 0 , 

(Xi,X2,X3) = < 

[ {+X3,b2,a2) if xi < 0 , 

in terms of constants Oj = ±1 and bi G R. By results in 0 ,y, 


( 8 ) 


18 l | , a system is locally approximated 


by dS]) when it satishes the conditions in ([2]). The brief outline of the dynamics that follows is not 
essential to the following sections, but we give it for completeness. 

The local flow ‘folds’ towards or away from the switching surface xi = 0, along the line X 2 = 
xi = 0 on one side of the surface, and along the line X 3 = xi = 0 on the other. Hence the point 
where these lines cross is called the ‘fwo-fold’. As a result, the surface xi = 0 is attractive in 
X 2 ,X 3 > 0 and repulsive in X 2 ,X 3 < 0 , while trajectories cross the surface transversely in X 2 X 3 < 0 . 
In the attractive and repulsive regions the flow slides along the surface xi = 0, and follows a vector 
held that is found by substituting ([ 8 |) into ([3]), and solving for A such that xi = 0. We will not 
discuss this sliding dynamics in detail, see for example [ 3 ] and references therein. 

The qualitative picture is then as shown in hgure [TJ The precise form of the local dynamics 
depends on whether the how curves towards or away from the discontinuity, determined by ai and 
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02, and also depends crucially on the quantity 6162, which quantifies the jump in the angle of the 
flow across the discontinuity. An accounting of the many classes of dynamics that arise from these 
simple conditions is given in j^, we give only the pertinent details here. 

The three main ‘flavours’ of two-fold are: the visible two-fold for oi = 02 = — 1 , the invisible 
two-fold for oi = 02 = 1, and the mixed two-fold for 0102 = —1; an example of each is shown in 
figure [T] (i,ii,iii) respectively. The terms visible or invisible indicate that the flow is curving away 
from or towards the discontinuity surface, respectively. In the cases depicted, there exist one or 
more canard trajectories, passing from the attractive sliding region to the repelling sliding region. 
This passage occurs in finite time (since the vector field obtained by substituting ([8]) into ([ 3 ]) is 
non-vanishing everywhere locally). The flow is unique in forward time everywhere except in the 
repelling sliding region, where it is set-valued because trajectories may slide along xi = 0, but 
may also be ejected into > 0 or xi < 0 at any point. This means that the flow may evolve 
deterministically until it arrives at the singularity by means of a canard, at which point it becomes 
set-valued, so we say that determinacy breaking occurs at the singularity whenever canards are 
exhibited. This occurs in the invisible case when 61,62 < 0 and 6162 > 1 , in the visible case when 
61 < 0 or 62 < 0 or 6162 < 1, and finally in the mixed case when 61 < 0 < 62 and 6162 < —1 or 
when 61 -|- 62 < 0 and 61 — 62 < — 2 . (The particular cases shown in figure [U are: (i) ai = 02 = —1 
with 61 < 0 or 62 < 0 or 6162 < 1; (ii) oi = 02 = 1 with 61,62 < 0 and 6162 > 1; (hi) 0102 = —1 
with 61 < 0 < 62 and 6162 < —1 or with 61 -|- 62 < 0 and 61 — 62 < —2.) 

This brief review follows the conventional picture, obtained by substituting ([8]) into ([ 3 ]). We 
shall only make the more general substitution of ([8]) into ([7]) when necessary, and will consider g 
to be a small perturbation, which will therefore have only a small effect on the dynamics outlined 
above; this can be studied more closely but is not our main concern here. Our interest now lies in 
what happens when we smooth out this system. 


III. REGULARIZING THE PIECEWISE SMOOTH SYSTEM 

We shall first outline the basics of regularization in a general form. Let x = {xi,X2,X3) and 
f = (/i) /2, /s)- We regularize the vector field in d?]) by replacing A with a smooth sigmoid function 

{ sign(xi) if Ixil > e 

> -h 0 (e) , > 0 for |xi| < e , ( 9 ) 

[-1,-M] if |xi| < e I 


for some small positive parameter e. In the Sotomayor-Teixeira regularization it is usually assumed 
more strongly that 4 >e{xi) = sign(u) for \xi\ > e, without the 0(e) term, but this not crucial to 
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our results. It is useful to introduce a fast variable u = xije. Noting by ([ 9 ]) that (peieu) = 
we obtain 


{eu,X 2 ,X 3 ) = {fi{eu,X 2 ,X 3 -(pi{u), f 2 {su,X 2 ,X 3 -(pi{u)), f3{eu,X2,X3]4>i{u))) 


(10) 


This is the regularization of ([ 7 ]). The region |xi| < e, which collapses to the discontinuity surface 
xi = 0 in the limit e —)• 0 , has been rescaled into the region |n| < 1 . In the following we restrict 
our attention to the dynamics in the region |u| < 1. 

Some basic properties of the system (uni) then follow using standard concepts of geometric 
singular perturbation theory (more detail of which can be found in j^, IT|, or various recent works 


in slow-fast dynamics such as 


22l | which will be relevant later). 


Expressing (IlOp system with respect to a fast time r = t/e, denoting the corresponding time 
derivative with a prime, gives the fast subsystem 

(u',X2,X3) = (/i(eu,X2,X3;0i(n)), ef 2 {eu,X 2 ,X 3 -,(j)i{u)), e/3(en, X2, X3; (/>i(u))) , (11) 

Setting e = 0 gives the critical fast subsystem (sometimes called the layer problem) 


(n',X2,X3) = (/i(0,X2,X3;(()i(n)), 0, 0) . 


( 12 ) 


This is a one dimensional subsystem, with a set of (x 2 , X 3 )-parameterized equilibria inhabiting a 
surface 


= {(u,X 2 ,X 3 ) G : fi{0,X2,X3;(pi{u)) = 0 ,\u\ < 1 } . 


( 13 ) 


This is an invariant manifold of (I 12 p wherever is normally hyperbolic, thus excluding the set 
of points where hyperbolicity fails, defined as 
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L = \ {u,X 2 ,X 3 ) € M : — /i(0,X2,X3;(/>i(u)) = 0 


( 14 ) 


Returning to (IlOp and setting e = 0 gives the slow critical subsystem (sometimes the called reduced 
problem) 


(0,X2,a:3) = (/i(O,X2,X3;0i(m)), /2(0, X 2 , X 3 ; (/>i(u)), f3{0,X2,X3]4>i{u))) 


( 15 ) 


which defines dynamics inside the critical manifold jC, on the original timescale. 

In our present context we will call Ad'^ the sliding manifold Ad'^, and refer to the dynamics in 
(fTHI) on Ad"^ as sliding dynamics, and to its solutions as sliding modes; this is due to their conjugacy 
to sliding modes in the piecewise smooth system as proven in [l 4 . . In Q it is shown that 
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the slow dynamics of (fTOj) on jC, is conjugate to the Filippov sliding dynamics of the piecewise 

smooth system ([ 3 j), and in [l^ this result is extended to include the generalization (j 71 ) . 

We now apply these ideas to the normal form of the two-fold singularity. The final step in our 
analysis will be to transform (IlOp into the normal form of a folded singularity, but we shall find 
that this is only possible if we regularize ([ 7 j) with g 7^ 0. 

IV. THE UNPERTURBED SYSTEM 

In this section we show that the Sotomayor-Teixeira regularization of the two-fold singularity 
has a sliding manifold as defined in section imi but that the non-hyperbolic line £ lies along 
the fast direction (the tt-axis), and hence in a degenerate position with respect to the flow. We 
show that this degeneracy is broken by perturbations of the form ([ 7 j) in section |Vl 

We first express the two-fold normal form as a convex combination by substituting ([8]) into ([ 3 ]), 
giving 

{xi,X2,X2,) = ai, 61)-h ^-^(0:3,62,02) ( 16 ) 

:= (/i(xi,X2,X3; A), f2{xi,X2,X3; X), /3(2:i, X2, X3; A)) . 

We then regularize this by replacing A i—>■ (p£{xi) for small e, constituting a Sotomayor-Teixeira 
regularization of and let u = xijs to obtain the two timescale system 

{eu,X2,X3) = - - -1(-X2,ai,6i) H- - - -{x3,b2,a2) (17) 

= {fi{eu,X2,X3;(pi{u)), /2(ew, 2:2, X3; (/>i(u)), /3(eu, X2, X3; (/>i(n))) . 

By (| 13 p . the sliding manifold is given by 

Ad'^ = I {u,X2,X3) G : (j)i{u) = — -|m| < 1 1 . ( 18 ) 

I X3 + X2 } 

The condition |i?i>i(ri)| < 1 by Q implies that Ad'^ exists only for 0 < X2X3 or for X2 = X3 = 0 . By 
expressing Ad'^ implicitly as the zero contour of the smooth function (X3 -|- X2)</>i(n) -|- (x2 — X3), 
we see that it is a smooth surface which twists over near X2 = X3 = 0 . It consists of two normally 
hyperbolic branches, one attractive in X2,X3 > 0 since dfi/du = —(X3 + X2)(p'i{u)/2 < 0 (using the 
fact that 4>'i{u) > 0 by Q), and one repulsive in X2,X3 < 0 since dfi/du = —(X3 + X2)(t>'i{u)/2 > 0. 
The two branches are connected at X2 = X3 = 0 along the non-hyperbolic set, found from (I 14 p to 
be 


£ = {(n, X2, X3) G Ad*^ : X2 = X3 = 0 } 


( 19 ) 
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This line segment £, at which the attracting and repelling branches of A 4 '^ intersect, constitutes 
the regularization of the two-fold singularity (xi = X2 = X3 = 0 in ([I])), now existing for all 
|tt| < 1 at X2 = xa = 0 on A 4 ^. Figure [ 2 ] shows an example of the piecewise smooth system (i), its 
regularization showing and L in (ii), which is then rotated in to show C more clearly (hi). 


regularize 


rotate 



FIG. 2: Regularizing the unperturbed system (fT6|l . for the example of an invisible two-fold, (i) The flow directions 
outside xi = 0 create an attracting sliding region in X 2 ,X 3 > 0 and repelling sliding region in X 2 ,X 3 < 0. (ii) 
The regularization of xi = 0, replacing xi with a fast variable u = xi/e, where the sliding regions create a 
critical manifold (shaded), hyperbolic except along the vertical line C, which aligns with the fast (double 
arrowed) u dynamics, (iii) The dynamics in the manifold is best viewed along the U 3 axis of rotated coordinates 

U2 = X2+ X3, U3 = X2- X3. 


We then have the main result of this section. 

Proposition 1. In the regularization of the normal form two-fold singularity the non- 

hyperbolic set L of the sliding manifold Ai^ lies everywhere tangent to the coordinate axis of the 
fast variable. 

Proof. The non-hyperbolic set C forms a line with tangent vector ec = ( 1 , 0 , 0 ) in the space of 
(m, X2, X3), which means it lies everywhere parallel to the fast u-coordinate axis of the two timescale 
system (diD. This is related to the fact that all derivatives of fi with respect to the fast variable 
u vanish along £, not only the first derivative dfi/du = — (xs -|- X2)4i'{u) = 0 which defines C 
as the set X2 = X3 = 0 (since </>'(«) yf 0 for |m| < 1 by ([ 9 ])), but also all higher derivatives 
fx/dvl = — (x3 -|- X2)(j)^{u) = 0 for any r > 1 . Thus this constitutes an infinite codimension 
degeneracy. 

□ 

In the literature on smooth two timescale systems, the connection of attracting and repelling 
branches of a slow invariant manifolds has been well studied, leading to a generic canonical form 
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and requisite non-degeneracy conditions as described in 2 ^. In the present notation, the non¬ 
degeneracy of along £ requires the conditions 


/i = 0 


^ - 0 
du — ^ ^ 


9fi ^ 0 A 0 

dx2,3 A U U 


( 20 ) 


the first three of which are satisfied on L as given by () 19 p , while the fourth is violated everywhere 
on £. Therefore the degeneracy of £ prevents us relating m to the canonical form of such 
singularities in slow-fast systems. 

In light of this problem, the papers [I 2 l . 


],Q,Q take a different approach to handling the system 
m- While the degeneracy above is not remarked on explicitly, the difficulties that arise from it 
are, and are tackled by re-scaling the local variables to prove that canards persist for perturbations 
within the Sotomayor-Teixeira regularization. Here we will instead permit perturbations that 
constitute a more general regularization, and in doing so we are able to obtain the canonical form 
satisfying (1201) . 


V. THE PERTURBED SYSTEM 


We will now show that a certain perturbation breaks the degeneracy of the system in the 
previous section. 

To achieve this, first observe that adding constant terms or functions of the coordinates 
(xi,X2,X3) to (fT 71 ) would only move the set £ in the {x2,x^) plane, not remove its degeneracy, 
easily seen since the derivatives 


^/l(0,X2,X3,(()i(u)) = /i(0-^2-^3,-U)^/i( 0,X2,X3, for r>0 

would still vanish on £, where /i( 0 , X2, X3, -|- 1 ) — /i( 0 , X2, X3, — 1 ) = 0 . The only recourse to break 
the degeneracy, specifically to give 7^ 0, is therefore to add terms nonlinear in 4 >i to (fTTll . and 
hence terms nonlinear in A to (I 16 p . Anything we add to the function fi in (I 16 p must still give ([8]), 
so it must vanish outside the switching surface xi = 0 , i.e. be a perturbation in the form ([ 7 ]). We 
shall show that perturbing xi with a small term proportional to — 1 is sufficient for structural 
stability. Perturbing X2 or x^ is neither necessary nor sufficient, therefore we leave them unaltered. 
The perturbed system we consider, applying ([ 7 ]) to the normal form ([8]) with g = (a, 0 , 0 ), is 


{xi,X 2 ,X 3 ) 


(-X2, ai, 61) (X3,62, 0-2) + (1 - A^)(a, 0, 0) 

{fl{xi,X2,X3;X)j2{xi,X2,X3;X)j3{xi,X2,X3;X)) , 


( 21 ) 
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where a is a constant. The regularization, replacing A for small e and letting u = xi/e, 

becomes 


{eU,X2,X3) 


^ (-X2, ai,bi) + - —(X3,62, 02) + a(l - 0?(u)) (1,0,0) (22) 

(/i(a;i,X2,X3; A), /2(xi, X2, X3; A), /3(xi, X2, X3; A)) , 


which is an a-perturbation of (I 17 h . Our main result is then: 


Proposition 2 . The regularization I 122 \) of the normal form two-fold singularity using & with 
g = (0,0,0), can be transformed into the canonical form for the folded singularity 12 BJ . namely 


£xi = X2 + X3 + 0 (exi,ex3,xix3) 
X2 = PX3 + gxi + 0 (x|,XiX3) 

X3 = r + 0 (x3,xi) 


provided a / 0 for small e > 0 , where p, q, r, are real constants, and provided the conditions 
i(6i — 62) < 1 = ai = —0,2 or ^{bi — 62) > —1 = oi = —02 do not hold. 

It turns out that the case excluded by the conditions ±^(61 — b2) < 1 = ±ai = =F“2 is that in 
which there are no canards or faux-canards, i.e. no orbits of the sliding flow passing through the 
singularity. 

We shall prove the proposition by way of three lemmas, establishing first the non-degeneracy 
of C, second locating a singularity along C, and finally using these to derive new local coordinates 
in which becomes a simple parabolic surface. 

The sliding manifold, found by applying (fT^ to (l 2 ^ . is now the set 

= |(m,X 2 ,X 3 ) G : |u| < 1, — ^i^X3 - t^^^X2 + a(l - (^i(n)) = o| , (23) 

which is normally hyperbolic except on the set given by applying (| 14 ll to ( 1221 ) . 

^ ^ j, , , ^ 2 a-bX 3 -X 2 X 3 -bX 2 'l 

C = < (u,X2,X3) : f>i(u) = 2 -^- =--- > . ( 24 ) 

I X3 -I- X2 4 a J 

Solving the conditions in ( 1241 ) we can express £ in paramteric form as 

(u,X 2 ,X 3 ) = £(yi):=(u,a((f>i(u)-lf,-a((/>i(u)-l-lf) ■ ( 25 ) 


This gives us the first result as follows. 

Lemma 3 . The non-hyperbolic set C, is transverse to the fast direction of i22\) . 
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Proof. By differentiating (I 25 p with respect to u, we find that the curve £ has tangent vector 
e£ = (1 , 2 a{ 4 >i{u) — l) 4 >'i{u) , — 2 a{ 4 >i{u) + l)(/)']^(u)), which for all \ 4 >i{u)\ < 1 is transverse to the 
coordinate axes provided a / 0. □ 

While the non-hyperbolic curve £ is now in a general position with respect to the fast variable, 
generically there may exist a new singularity along £, where the flow’s projection along the A- 
direction onto the nullcline /i = 0 is indeterminate. This is the so-called folded singularity, defined 
in the following lemma. 

Lemma 4 . For the values of the constants 01,02,61,62, given in Proposition there exists an 
isolated singularity of the flow along the non-hyperbolic set C, where the projection of the slow flow 
onto lies tangent to £. 


Proof. Let us consider the slow critical subsystem, obtained by letting e = 0 in (| 22 ] 1 . 

{ 0 ,X 2 ,x^) = ^^^^-^^{-X 2 ,ai,bi) + - — j^{x 3 ,b 2 ,a 2 ) + a{l-(j)l{us)){l, 0 , 0 ) ( 26 ) 

= (/l(O,X2,X3;0i(u)),/2(O,X2,X3;0i(u)),/3(O,X2,X3;())i(u))) . 


Since is the surface where /i = 0 , a solution of (I 26 p that remains on for an in¬ 
terval of time satisfies /i = 0 . We can find ii on A 4 ^ using the chain rule, writing /i = 
(u,X2,X3) ■ {d/du,d/dx2,d/dx3) fi = + {f 2 , h) ■ d{il,x3) = which rearranges to o = 

~(/2,/3) • J^xl,x3) / (^&)' Thus u is indeterminate on A 4 ‘^ at points where the numerator and 
denominator of this vanish, or in full, where 


OU d{X2,X3) 


( 27 ) 


These three conditions define an isolated singularity on £ C Ai^. Denoting the value of fl at the 
singularity as fis, and solving which constitutes finding a point on £ given by (IMIl such that 

l + (pl{Us) 


0 — if 2 s, / 3 s) • - 


_ ( ai+b2 1 ai-b2 

2 w 2 


2 

4 >i{us) , 


bi+a2 I fci—Q2 

2 w 2 




l + 4ll{Us) 1 —(/>l(us) 
2 ’ 2 


( 28 ) 


we find that the folded singularity lies at («,^2,^3) = {us,X2s,X3s) where 

ai+a2 £ I 4aia2 

4 >l{Us)= > X 2 s = a{ 4 >i{Us) - , XSs = -a{ 4 >l{Us) + if . ( 29 ) 

^ + bl-b 2 

Noting that ai and 02 in the normal form just take values ± 1 , and recalling that 4 >i{us) is monotonic 
on |rt| < 1 by Q, we have: 
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• in the case oi = 02 = 1, we have (t>i{us) = =t y 1 + ; implying that there exists 

a unique solution Us G (—1, +1) for any bi and 62 (the positive root for hi >62, the negative 
root for 61 < 62); 

• in the case ai = 02 = —1, we have (t)i{us) = 573^ =t y^l + , implying that there exists 

a uniquesolution Ug G (—l,+l)for any 61 and 62 (the positive root for hi < 62, the negative 
root for 61 > 62); 

• in the case ai = —02 = 1, we have (/>i(us) = implying that there exist two 

solutions Us G (— 1 , + 1 ) for 61 — 62 > 2 , and no points otherwise. 

• in the case oi = —02 = —1, we have 4 >i{us) = implying that there exist two 

solution Us G (— 1 , + 1 ) for 61 — 62 < — 2 , and no points otherwise. 

□ 


This lemma establishes the existence of at least one unique folded singularity on C in the cases 
listed in Proposition [ 2 j In the cases where Us is unique we proceed directly to the steps that follow 
below. In the cases where Ug can take two values we can proceed with the following analysis about 
each value, and will obtain different constants in the final local expression, i.e. a different folded 
singularity corresponding to each Ug. In the cases when Ug does not exist, no equivalence can 
be formed; these are the cases when the the two-fold’s sliding portrait is of focal type (see j^), 
and there exists no canards since orbits wind around the two-fold but never enter or leave it. So 
excluding those cases oi = —02 = I with 61 — 62 < 2 and ai = —02 = — I with 61 — 62 > —2, we 
proceed with the final step in proving proposition [2l 

Lemma 5 . Coordinates can he defined in which the folded singularity of ( 2 ^ lies at the origin, 
and £ lies along a coordinate axis corresponding to a slow variable. 

Proof. Taking a valid solution of Ug from (| 29 l) for |u| < I, a translation puts the singularity at the 
origin of the new coordinates 


yi = (j)l{u) - (fiiug) , y2 = X2-X2g, y 3 = X 3 -XSg. 


( 30 ) 


Then fi becomes 


l + 4>l{Ug) l-fil{Ug) 

h = --2/2 H--2/3 


2/3 + 2/2 


+ ayi yi 


2 


2 


2 


( 31 ) 
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found by using (I 28 |) - (f 2 ^ to ensure that terms involving X2s and x^s vanish. To find coordinates in 
which L lies along a coordinate axis, from (12511 we can obtain the yi-parameterized expression for 


(yi,y 2 ,y 3 ) = ^(yi) := (yi,-ayi (2 - - yi), -ayi{2 + 24 >i{us) + yi)) , 


and re-arrange this to take ys as a parameter, expressing C as (yi,y2) = {yiL{yz)iy2L{yz))-, where 
yiL(yz) \ I -1 - -b 

y2L{yi) ) \ -ys - 4 a(-i - ^i{us) + {i + 4 >i{us)y -ys/a) 

The derivatives of these functions are needed to evaluate the vector field components below, these 
are 



y'ldys) 


-l/2a 

I + (j)i{us) + yidys) ’ 


ydiys) 


1 - (j)i{us) -yiL(y 3 ) 

1 -b (j)i{us) +yiL(y 3 ) ' 


We can then rectify £ to lie along some axis by defining new coordinates 


( 33 ) 


zi = yi- yidys ), ^2 = y2 - y 2 L(y 3 ), ^^3 = y 3 • ( 34 ) 


The original vector field components can then be written as 

^ _ ^,2 _ ,^, J 2 , 

f2 = f2s + (Zl + yidzs)) = / 2 s + 0 (Zi, Zs) , (35) 

fs = /Ss + {zi + yidzs)) ^ = hs + 0 {zi,Z3) . 

With a little algebra we find that 


ezi 


eyi - ^my'iLiys) 

l + (j)i{us) ( efss 

2 Va(l + ())i(us ))2 


Z2dl{Us 


- azldius) + 0 [sZ3,eZi,Z2Z3,Z2Zi,zf) 


A small shift Z2 = Z2 — a{i+(t>i{usd(t>' after some lengthy but straightforward algebra, 

using the relations in (j 28 p and () 35 p to show that any terms not propotional to zi or Z3 vanish. 


Z2 = f2- zsydizs) 

P 

= qzi + -Z3 -b 0 (Z3,ZIZ3) 
a 


df 2 s dfssl - (f)i{us) _ 2 / 3 s + g(l-b (/>i(us))^ 

d 4 >i d(j)i 1 + (pi{us) ' ^ 2(1 + (/)i(us ))2 


where 
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The last thing to do is just scaling. Collecting everything together so far we have 


£Zi = (diZ2 - azf)(/)[(Us) + 0 (£Z,£Z3 ,ZiZ3) 

^2 = ^Z3 + qzi + 0 (zl, Z1Z3) 

Z3 = hs + 0{z3,Zi) 

where di = —^(1 + ct)i{us))- Defining new variables xi = ■s/\a\<j)'i{us)zi, X2 = —sign{a)di(j)[{us)z2, 
X3 = —sign(a)2;3, and t = —sign(a)t, gives 


£Xi = X2 + xf + 0 {£Xi,£X3,XiX3) 

X 2 = PX 3 + qXl + 0 [xl,XiX 3 ) ( 36 ) 

X3 = f + 0(x3,Xi) 


where 


( 37 ) 


f = h,, P = {k- + A, - 29Vlak'iK)) , 

^ = vnlife) (*'^‘*"‘> + -1) wf) ■ 

Omitting the tildes, this is the result in the lemma and in Proposition [ 5 l clearly valid only for 
a 7^ 0. 


□ 


Remarks on the singularity 

A glance at the papers 
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2 , 'll ] reveals what a charismatic singularity lies hidden in the 


dynamics of the two-fold, waiting to be released when the piecewise smooth system is perturbed 
by simulations that smooth, regularize, or otherwise approximate the discontinuity. As for the 
two-fold itself in section m a detailed description is beyond our interest here and can be pursued 
in future work, but as a guide we shall briefly gather together the main points from the literature. 




For sufficiently small e > 0 , by standard results of geometric singular perturbation theory 
there exist invariant manifolds in the neighbourhood of A 4 '^/£, on which the dynamics is 
topologically equivalent to the sliding dynamics found above. Trajectories that pass close to the 
singularity, or more precisely, close to the folded singularity on the set C, may persist in following 
the manifold from its stable to unstable branches, while other nearby trajectories will veer 
wildly away, their fate sensitive to initial conditions and proximity to primary canard orbits (those 
which persist along both branches of Ai^ throughout the local region). 

Figure [ 5 ] shows an example of the perturbed system and its regularization for each flavour of 
two-fold in (i) (corresponding to those in figure [T]), followed by their regularization (ii), and a 
rotation (iii) to show the phase portrait around the set L more clearly (similar to figure [ 2 ]). In the 
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most extreme case, the folded node, the original phase portrait contains infinitely many intersecting 
trajectories traversing the singularity, while the perturbed system splits these into distinguishable 
orbits, a finite number of which asymptote to the attracting and repelling branches of the critical 
manifold. 



FIG. 3: blowing up the perturbed (a 7 ^ 0) system, for examples of each flavour of two-fold. Labelling as in 
figure[2] Note in the regularization (ii) that £ is now a curve. Rotating around the u axis in (iii) we can see the 
attracting branch (upper right segment) and repelling branch (lower left segment) of the sliding manifold A4^ 
(shaded), connected by £. The folded singularity (f.sing.) appears along £, two in the case of mixed visibility, 
recognised as having a phase portrait that resembles a saddle or node if we reverse time in the repelling branch 
of In (iv) we sketch the corresponding phase portraits in the slow-fast system ()36p . 

Like the different kinds of two-fold, there are different classes of folded singularity, and their 
classification depends on the slow dynamics inside . From the expressions (I 36 ]l - (l 371 l we see that 
the class therefore depends not only on the constants oi, 02, 62, of the original piecewise smooth 

system, but also on the ‘hidden’ parameter a. 

The classification scheme is fairly simple, and can be used to verify the dynamics on seen 
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in figure [3l The projection of the system (I 36 p onto A 4 ^, found by differentiating the condition 
0 = X2 + xf with respect to time to give 0 = 6x3 + cxi + 2xixi + 0 (x|, X1X3), is 




+ 0 (x|,XiX3) 


A classification then follows by neglecting the singular prefactor l/ 2 xi and considering whether 
the phase portrait is that of a focus, a node, or a saddle. This is determined by the 2 x 2 matrix 
Jacobian, which has trace c, determinant 2 db, and eigenvalues ^(c± \/(? — 8db). This will not 
be the true system’s phase portrait because the time-scaling from the l/2xi factor is positive 
in the attractive branch of , negative (time-reversing) in the repulsive branch, and divergent 
at the singularity (turning infinite time convergence to the singularity into finite time passage 
through the singularity). The effect of this is to ‘fold’ together attracting and repelling pairs of 
each equilibrium type, so each equilibrium becomes a ‘folded-equilibrium’, forming a continuous 
bridge between branches of . 

As a result the flow on is a folded-saddle if db < 0, a folded-node if 0 < 8db < (?, and a 
folded-focus if (? < 8db. Canard cases occur for c > 0 and faux canard for c < 0 . In figure[ 3 ]we show 
the result of regularizing the piecewise smooth system for an example of each type of two-fold that 
exhibits determinacy-breaking (those from figured]). In the visible two-fold the singularity becomes 
a folded-saddle, in the invisible case it becomes a folded-node, while the mixed case becomes a pair 
consisting of one folded-saddle and one folded-node. 

One may ask why certain cases were excluded by the proposition above. The excluded cases 
were those in which no canards exist in the slow-fast system. Canards occur when transversal 
intersections exist between the attracting and repelling branches of the slow manifolds. If no 
such intersections exist, the critical system possesses no folded singularities and hence is excluded 
from Proposition [21 Hence the omission of these cases is consistent, and a posteori it is obviously 
necessary, in the equivalence sought in the proposition. 

With our main result proven, we conclude with two sections relevent to the study above, which 
help elucidate certain ideas that have arisen lately in the study of piecewise smooth dynamical 
systems, which are particularly relevant to the study of singularities like the two-fold. 
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VI. “DISCONTINUITY BLOW-UP”, AN APPROACH TO NONSMOOTH SYSTEMS 

WITHOUT REGULARIZATION 

Many of the relations in the analysis above had to written implicitly in terms of the regular¬ 
ization function 4 >i{u), having introduced A = (peixi) = (piiu). We could have preceded instead 
by studying the dynamics of A directly, omitting reference to a smoothing function altogether. 
This alternative approach to nonsmooth systems was discussed in Q], and makes the results above 
somewhat more concise and explicit, but a couple of propositions are required to show that stan¬ 
dard concepts from geometric singular perturbation theory can be applied on {X,X2,x^) in place 
of {u,X2,X3). Having followed the conventional route above, we provide the basic results needed 
for this alternative route here. 

Firstly we require a dynamical system on A. 

Proposition 6 . The dynamics of A is given by 

e(A, e)A = xi = f(x, A) • Vxi ( 38 ) 

such that e(A, e) <C 1 , where e denotes a continuous positive function and e a small parameter, 
with 0 < e < e* <C 1 and A G terms of constants e* and that satisfy ^ ±1 as 

e* ^ 0. 

Proof. Consider regularizing the vector field ([ 7 ]) by replacing A with a differentiable sigmoid func¬ 
tion 4 >s{xi) as defined in ([ 9 ]). We shall use the relation A = 4 >e{xi) to derive a dynamical system 
on A. Differentiating A = with respect to t gives 

A = —())(.(xi) for \xi\ < e . ( 39 ) 

Considering a variable u = xije we see that, according to ([ 9 ]), the function <^£(xi) = <^i(n) and 
its derivative £())(.(xi) = (f'liu) are smooth with respect to u in the limit e ^ 0 . Moreover 4 >i{u) is 
strictly positive because <^i(u) is strictly increasing, and (/>'^(u/e) only becomes small (or vanishing) 
for |xi|/e > 1 . So the quantity s/cp'i is small and nonzero for |xi|/e < 1 , and using it we dehne 
a fast timescale r = t(/>(.(xi)/e. Since 4 >i{u) is differentiable and monotonic for |a| < 1, it has an 
inverse ipiX) such that = u, and we can define a function 

e{X,e) := s/(f[{fj{X)) , for |A| < 1 . ( 40 ) 

That this quantity is small is shown as follows: the function (pi{u) varies differentiably over interval 
on which its extremal values are = ±1, therefore there exists a point tt* where (j)'[u^) = 
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= Ij by continuity since (^>'^(±1) = 0, there exist two points uf where 4)'{uf) = =be 
for 0 < e < 1 , and moreover an interval u~ < u < uf such that 4 >iiu) > e. Fix some e* such that 
0 < e* ^ 1, then s/(j)i{u) < 1 for < u < , and 

lim e(A, e) = 0 , 

£^0 

so that e(A,e) <C 1 for e <S e* and u G 

By (f 39 ]l we therefore have the dynamical equation e(A,e)A = xi = f(x, A) • Vxi in the proposi¬ 
tion. □ 

This proposition identifies A as a fast variable inside A G (— 1 , -|- 1 ) (more strictly for A G 
where and e* is arbitrarily small but nonzero). When A is set-valued on n = 0 with 

e = 0, this equation determines its variation on the timescale r which is instantaneous relative 
to the timescale t. In the piecewise smooth dynamics literature this is sometimes referred to as 
providing the blow-up of the discontinuity surface xi = 0 into an interval A G [— 1 , -|- 1 ]. 

The system obtained by applying (l 38 l) to ([ 8 ]), using ([ 7 |) on the interval A G then 

{eX,X2,X3) = {fi{xi,X2,X3;X), /2(xi, X2, X3; A), h{xi,X2,X3 ;X)) . ( 41 ) 

Where possible we can omit the arguments of e = e(A,e) without confusion. While standard 
geometrical singular perturbation theory does not apply to (|^T]1 because e is a function, we can 
show easily that this leads to the same critical manifold geometry as the conventional approach 
outlined in section HU 

Proposition 7 . The system 0 has equivalent slow-fast dynamics to the system m on the 
discontinuity set xi = 0 in the critical limit e = 0. 

Proof. Rescaling time in (l 4 T]l to r = t/e, then setting e = 0 and xi = 0 , gives the fast critical 
subsystem 

(A',X2,X3) = (/i( 0 ,X 2 ,X 3 ; A), 0 , 0 ) . ( 42 ) 

The equilibria of this one-dimensional system form the manifold 

= {(A,X2,X3) G [- 1 ,-M] X : /i( 0 ,X2,X3; A) = O} , 

which is equivalent to the manifold () 13 p . This is an invariant manifold of the system () 4 ip in the 
e = 0 limit everywhere that is normally hyperbolic, that is excepting the set 

C = |(A,X2,X3) G : ^/i( 0 ,X 2 ,X 3 ; A) = 0 
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Since l^/i = </.;(«) §xfi and (p'i{u) / 0 for |it| < 1, this definition of C is equivalent to (fHl) . Setting 
e = 0 and xi = 0 in (HIT) gives the slow critical subsystem 

( 0 ,X 2 ,X 3 ) = (/i( 0 ,X 2 ,a: 3 ; A), /2(0, X2, X3; A), /3(0, X2, X3; A)) , ( 43 ) 

which defines dynamics in the critical limit e = 0 on which is exactly as given by (USD. 

□ 

In Q] an extension to the Filippov approach to piecewise smooth dynamics is therefore proposed, 
introducing directly a dummy timescale t = t/e. Denoting the derivative with respect to this fast 
time r we have simply 

A^ = xi = f(x, A) • Vxi , 


and in full 


{X',X2,X3) = {fi{xi,X2,X3;\), /2(X1,X2,X3; A), /3(X1,X2,X3; A)) . ( 44 ) 

This permits study of the slow-fast critical dynamics without reference to regularization functions, 
providing a more direct analysis of sliding dynamics. The system ( 144 ^ directly permits us to 
consider the discontinuity surface xi = 0 , and fixes the value that A takes on the interval [— 1 , + 1 ] 
inside the manifold (when it exists) for which the set xi = 0 is invariant, and specifying the 
variation (x2,X3) on that manifold. In regions of xi = 0 where does not exist (where /i 7^ 0 
for xi = 0 and A € [— 1 ,+!]), the fast subsystem of (l 44 l) conveys the flow across the discontinuity 
from one region xi ^ 0 to the other. This analysis via the dummy system () 44 ll can be applied 
to reformulate the previous sections without regularization, yielding equivalent results, the main 
difference being that <?!>i(u) and (i)'i{u) can be replaced everywhere by A and I respectively. 

VII. NONLINEAR SWITCHING TERMS, AND THE AMBIGUITY OF 

REGULARIZATION 

The vector held in section IVl demonstrates that qualitatively different geometry can be obtained 
if we consider perturbations outside the Sotomayor-Teixeira regularization. The dynamics that 
results is rather complex, however, so we shall briehy show an example of two non-equivalent 
smooth systems with the same piecewise smooth limit &■ This highlights that there must be 
some ambiguity in how to regularize the discontinuous system into a smooth system. We show 
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that the nonlinear dependence on A in ([ 7 ]) provides a means to disctinguish between the different 
possible regularizations, and study them in the piecewise smooth limit. 

Consider two smooth system given by 


{x,y) = {(pe{x),l) , 

( 45 ) 

{x,y) = {(pe{x),l -2(pl{x)) , 

( 46 ) 


where e is a small positive parameter and <^e(x) is a sigmoid function as defined in ([ 9 ]). We will 
refer to () 45 p and ( 1461 ) as the linear and nonlinear systems, respectively. 

A common approach taken to studying sigmoid systems like (| 45 p or (| 46 p in applications, par¬ 
ticularly when the systems above represent empirical models, is to study the system obtained by 
replacing (psix) with sign(x) in the limit £ = 0 . The result is a piecewise smooth system given for 
X / 0 by 

(x,?/) = (sign(x),-l) , ( 47 ) 

for both (I 45 p and (jUj). The system is set-valued at x = 0 , and the question is then how to solve 
the differential inclusion at x = 0 . In the Filippov method we replace sign(x) with a switching 
parameter A, such that A = sign(x) for x 7^ 0 and A G [— 1 , -|- 1 ] for x = 0 , and write 

(x,7/) = (A,-l) , A G [- 1 , 4 - 1 ]. ( 48 ) 

This is simply the system obtained by applying the convex combination ([ 3 ]) to (HTl) . A figure HI 
illustrates, the two vector fields point towards x = 0 so the entire discontinuity line must be 
a sliding region. The Sotomayor-Teixeira regularization smooths this by replacing A again with 
another sigmoid function 4 >£{x) (obeying the same properties as in ([ 9 ])) 

{x,y) = {(p^ix),-!) . 

This recovers the system (I 45 p . but not () 46 p . The dynamics of the two is crucially different. The 
linear system ( 145 p has a critical manifold, given by n = 0 if we assume <(>1(0) = 0 . The dynamics 
on the manifold manifold is given by {u,y) = ( 0 ,- 1 ). To see this, introducing a fast coordinate 
tt = ex, we have 


{£u,y) = {(pi{u),-l) 

and apply standard geometrical singular perturbation [s, llj in the limit e — 5 - 0 . For small positive 
£ the dynamics near u = 0 should be a small perturbation of {u,y) = ( 0 , — 1 ). This is equivalent 
to the Filippov sliding dynamics on x = 0 in the piecewise smooth system, (x,y) = ( 0 , — 1 ). 





21 


We failed to obtain (H6l) under regularization because we omitted what happens to the 
term in the piecewise smooth limit (as e ^ 0). If we instead apply ([7]) to (jl7|l with g = (0, 2), we 
obtain 


(x.y) = (A,l-2A2). 


(49) 


When we regularize this we obtain 

{x.y) = {(l)e{x),l- 2(/)2(x)) , 


and thus we have regained the nonlinear system (I46p under regularization. The nonlinear system 
again has a critical manifold, given by u = 0 if we assume <^i(0) = 0, but now the dynamics on 
the manifold manifold is given by (u, y) = (0,+1). This is seen by introducing a fast coordinate 
u = £x, giving 

i£u,y) = i4>i{u)A - 20?(u)) 


and applying standard geometrical singular perturbation theory in the limit e ^ 0. For small 
positive £ the dynamics near n = 0 should be a small perturbation of {u,y) = (0,+1), and not in 
any way close to that of the linear system (which points in the opposite direction), as illustrated in 
figure m This dynamics is also not equivalent to the Filippov sliding dynamics of (I47h via (I48h . but 
is instead equivalent to the nonlinear or ‘hidden’ sliding dynamics of (1471) via (1491) . as described in 





FIG. 4: Two approaches to the piecewise smooth system dUl): the linear combination (j48p and nonlinear 
combination (jH]), with their regularizations (j43|) and (|46]). 


VIII. CLOSING REMARKS 


Regularizing or smoothing a discontinuity of course raises issues of uniqueness, namely that 
infinitely many qualitatively different smooth systems can have the same piecewise smooth limit. 
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The idea of nonlinear switching terms — nonlinear dependence on the parameter A in ([7]) — is 
that they provide a way of restoring uniqueness by distinguishing the limits of different smooth 
systems. We have shown here that the smooth system obtained from a two-fold singularity is 
structurally unstable if it depends only linearly on A, but that a small perturbation, by terms that 
are nonlinear in A, restores structural stability and allows transformation into the general local 
singularity expected in a smooth system, namely the folded singularity. 

For such a simple system, even taking its piecewise linear local normal form, the two-fold 
exhibits intricate and varied dynamics. Just how intricate becomes even more clear as we attempt 
to regularize the discontinuity, and study how the two-fold related to slow-fast dynamics of smooth 
systems. As well as insight into the dynamics that is seen upon simulated such a system, this 
adds a new facet to the question of the structural stability of the two-fold, which has remained a 
stimulating question since [l^ . 

An in-depth description of the dynamics that ensues in the different cases of two-folds, and the 
smoothings subject to perturbations, would be lengthy, and deserves future study elsewhere. As 
a demonstration we conclude with three examples showing the complex oscillatory attractors that 
can be formed by two-fold singularities. We take 

(i) f+ = (-X2, fiCi -b ^X2 - 1, ^X2 - \x 2 X 3 - |), f” = (X3, \x2Xz - f, |x3 - 1 - Xi)' 

(ii) f+ = (-X2,1 -b Xi, -|), = (X3, 1 - fxi); 

(iii) f+ = (-X2 -b ^xi,xi - |,xi - 2), f" = (x3 -b -b 1 - xi); 

and use [7] with a = 1/5, and simulate them in figure [5j 

The numerical solutions apply Mathematica’s NDSolve to a regularized system, replacing A = 
sign xi by a sigmoid function tanh(xi/e) with e = 10“^. Further simulations not shown here verify 
that such dynamics persists with different monotonic smoothing functions, such as xxj\J -b xf, 
and for different values of e. The first row shows the simulation of a single trajectory for a time 
interval t = 1000 in (xi,X2,X3) space, as the flow attempts to switch between the the xi > 0 and 
xi < 0 flows, and the critical slow manifold flow. The piecewise smooth system is sketched in 
the second row, and the regularization in the third row, including the critical manifold . The 
sliding vector field has canard trajectories passing from the righthand attracting branch to the 
lefthand repelling branch via the folded singularity; at a visible two-fold only one canard exists, 
at an invisible two-fold every sliding trajectory is a canard, and at a mixed two-fold a region of 
trajectories are canards. 
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FIG. 5: Three examples of attractor organised around a two-fold singularity. Showing: (i) a simulated trajectory, 
(ii) a sketch of the piecewise smooth flow inside and outside xi=Q that gives rise to it, and (iii) the blow up on 

xi = 0 . 

The degeneracy in section IIVI gives some insight into why simulations of systems containing 
two-fold singularities with determinacy-breaking singularities exhibit highly sensitive dynamics, by 
relating it to a canard generating singularity in a singularly perturbed system. 

In the ongoing saga of the two-fold, the system (|2T|) now succeeds (fTBIl as our prototype for 
the local dynamics. The question of whether this constitutes a ‘normal form’ has issues both in 
the piecewise smooth and slow-fast settings, but it is clear that (1^ is structurally stable, and 
represents all classes of behaviour that occur both in the piecewise smooth system, and in its blow 
up to a slow-fast system. 
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